ELRADGEN: Monte Carlo generator for radiative 
events in elastic electron— proton scattering 

A.V. Afanasev 

Jefferson Lab, Newport News, VA 23606, USA 

I. Akushevich 

rn " 

§ : Duke University, Durham, NC 27708, USA 

<N 

Ojt)' A. Ilyichev 

=S ' 

<^ ; Center of Particle and High Energy Physics, 220040 Minsk, Belarus 



O 



B. Niczyporuk 
Jefferson Lab, Newport News, VA 23606, USA 



> 

O ■ Abstract 

We discuss the theoretical approach and practical algorithms for simulation of ra- 
1 diative events in elastic ep-scattering. A new Monte Carlo generator for real photon 

emission events in the process of elastic electron-proton scattering is presented. We 
"*^h , perform a few consistency checks and present numerical results. 

& 1 ; 1 Introduction 

Radiative effects accompany any processes with charged light particle scattering. These 
effects are exactly calculable in QED and they should be incorporated as an essential step 
in the data analysis. For elastic processes, which are the subject of the present report, our 
starting point for the radiative effect calculation is the technique proposed by Bardin and 
Shumeiko [I] and later developed for polarized particles as well as for specific elastic and 
exclusive processes in papers [21IH1SE]- For other approaches, please see Refs. [H10IH]- In 
some cases the radiative effects can be calculated analytically. However, often the analytical 
integration over the photonic phase space is complicated because of specifics of detector 
geometry or is not desirable because of experimental setup (for example, detector resolution 
and kinematic cuts). In these situations an approach based on Monte Carlo simulation of 
the radiative effects is the most adequate. Monte Carlo generator RADGEN [Oj for deep 
inelastic processes was developed on the basis of corresponding analytical code POLRAD 
10 and successfully applied by several experimental collaborations analyzing polarized and 
unpolarized observables in deep inelastic electron-nucleon scattering. 

The reason for careful implementation of the radiative correction procedure in elastic 
ep-scattering is twofold. First, elastic ep-processes are commonly used for detector cal- 
ibration, and JLab CLAS detector is one of the examples. Second, a new generation of 
experiments on elastic ep-scattering provide data with very high accuracy (see, for instance, 
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Refs. |11U12| ). and set new precision requirements to data analysis and, in particular, to 
the radiative correction procedure. In this report we present a new Monte Carlo generator 
for elastic unpolarized ep-scattering events taking into account QED radiative effects. The 
main requirement for construction of the generator is appearance of simulated events (Born 
and radiative) as well as all event distributions with respect to the kinematic variables in 
accordance with their contribution to the total cross section. 

The contribution to the total cross section from the Born process 



is described by the transfer momentum squared Q 2 = — (ki — k 2 ) 2 and the azimuthal angle 
(p. The cross section does not depend on the angle up to acceptance effects. To describe 
the phase space of the radiative process 



three new kinematic variable have to be considered. We choose them as the proton transfer 
momentum squared t = — (ki — k 2 — k) 2 , the inelasticity v = (p2 + k) 2 — M 2 and the azimuthal 
angle 4>k between the planes (q, k x ) and (k 1; k 2 ). This set of the variables defines the four- 
vectors of all final particles in any frame. Usually Lab system is used with OZ axis along 
the beam direction and plane OZX parallel to the ground level. 

Thus, the strategy for simulation of one event can be defined as follows. 

• For the fixed initial energy and Q 2 the non-radiative and radiative parts of the total 
cross section are calculated. Note, that together with Born subprocess, the contribu- 
tions of loop effects and soft photon radiation are also included into the non-radiative 
part. 

• The channel of scattering is simulated for this event in accordance with partial contri- 
butions of these two positive parts into the total cross section. 

• For any process is simulated uniformly from to 2n. 

• For the radiative event the kinematic variables t, v and <fik are simulated in accordance 
with their calculated distribution (see discussion below). 

• The four-momenta of all final particles in the required system form are calculated. 

• If the initial Q 2 has not a fixed value but is instead simulated according to some 
probability distribution (for example, the Born cross section) then the cross sections 
have to be stored for reweighting. The Q 2 distribution is simulated over the Born cross 
section, and realistic observed Q 2 distribution is calculated as sum of weights, they are 
ratios of the total and Born cross sections. 

We note that separation of the bremsstrahlung process into the radiative and nonradiative 
parts necessarily requires to introduce an additional separating parameter. Usually this 
quantity is associated with photon energy resolution in the detector. In our program the 
minimal inelasticity v m i n is chosen as this parameter. 
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Figure 1: Feynman graphs contributing to the non-radiated (a, b, c) and the radiated (d, e) 
parts of the cross sections in lepton-nucleus scattering. 

2 Non-radiated and radiated contributions 

The radiative events are generated in accordance with their contribution to the observed 
total cross section, that can be separated into non-radiative and radiative parts: 
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Here and later we define a = da/dQ 2 d<f). The first part includes the Born cross section o~q 
(Fig.l (a)), loop effects (Fig.l (b,c)) contributions as well as soft photon radiation, which is 
restricted by the inelasticity value v < v min : 
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The explicit expressions for o~o, ^m/, 8vr can be found in jHj. The dependence of the non- 
radiative part of the cross section on w min appears in the following way: 
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Here S = 2kxPi, v max = S — Q 2 — M 2 Q 2 /S, J-i(t) are the squared combinations of the 
electric and magnetic elastic form factors. The quantities 9i as well as Of are the analytical 
functions j3j of kinematic invariants for the radiative and Born subprocesses, respectively. 

The radiative part of the total cross section can be presented as an integral over the real 
photonic phase space in the following way: 
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where 
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Figure 2: The region of integration over v- and t- variables for JLab kinematics (Q 2 = 3 
GeV 2 , S = 7.5 GeV 2 ). 
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The radiative photonic phase space t- and v- variables are presented in Figure El 

Notice that in © as well as in (J7J) the structure functions Ti depend on only one inte- 
gration variable, namely, on t. Therefore to apply our program to different fits or models of 
nucleon form factors in (JHJ) and 0, the numerical integration over the variable t is used. On 
the other hand, to speed up the process of event generation, integration of the expressions 
0i in (JBJ) over the other photonic variables v and (fik has been performed analytically because 
the quantity 9i(t,v,(pk) appears as a result of the calculation of the matrix element squares 
with the real photon emission and do not depend on any model for elastic electron-nucleon 
scattering. 



3 Simulation of radiative events 

To simulate v and <pk in accordance with the analytically calculated distributions without any 
numerical integration over t, the first step is simulation of the variable t. Then taking into 
account the simulated value t, the distribution function for v can be constructed. Finally, 
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Figure 3: t- and v -histograms and the corresponding probability densities for the kinematics 
Q 2 = 3 GeV 2 and S = 7.5 GeV 2 . 

<p k is simulated according to the simulated values of t and v. If the event generation is 
performed successfully, the event distributions over these three variables has to follow the 
corresponding probability densities: 
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As an example, the simulated distributions of the kinematic variables t and v together 
with their probability densities for representative kinematics at CEBAF energies are com- 
pared in Figure 3. 

In summary, we presented the first results of the FORTRAN code ELRADGEN, which is 
a Monte Carlo generator of elastic unpolarized ep-scattering that includes electromagnetic 
radiative events. 
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